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We study the contribution of the strong interactions between the two pions in S- and P-waves to 
the weak B — > tt-kK decay amplitudes. The interference between these two waves is analyzed in the 
7T7T effective mass range ofthep(770)°and/o(980) resonances. We use a unitary irn and KK coupled- 
channel model to describe the 5-wave interactions and a Breit-Wigner function for the P-wave 
amplitude. The weak B-decay amplitudes, obtained from QCD factorization, are supplemented with 
charming penguin contributions in both waves. The four complex parameters of these long distance 
terms are determined by fitting the model to the BaBar and Belle data on B*' — > n + ir~ K ±,a 
branching fractions, CP asymmetries, nn effective mass and helicity- angle distributions. This set 
of data, and in particular the large direct CP-asymmetry for B —> p(770)° K decays, is well 
reproduced. The interplay of charming penguin amplitudes and the interference of S- and P-waves 
describes rather successfully the experimental S and A values of the CP-violating asymmetry for 
both B° -> /o(980)Jf§ and B° -> p(770)°K% decays. 



PACS numbers: 13.25.Hw, 13.75Lb 
I. INTRODUCTION 

Hadronic final-state interactions in three-body charm- 
less B meson decays supply good opportunities for CP 
violation searches. In order to interpret in the most 
reliable way decay observables, it is necessary to take 
into account strong interaction effects between the pro- 
duced meson pairs. These effects, at relative energies 
below and about I GcV, are clearly visible in the Dalitz 
plots of three-body weak decays such as B — > irirK and 
B — » KKK [1—6]. In particular, one observes a distinct 
surplus of events at these relatively small effective masses 
distributed near the edge of the Dalitz plot. Prominent 
maxima are seen in the effective 7T7t mass distribution, 
notably the /o(980) and p(770)° resonances. Several res- 
onances at higher energies are visible too, though their 
identification is less straightforward. The experimental 
effective mass distributions of three-body _B-decays are 
extracted from Dalitz plots. One cannot describe these 
data without including final state interactions. The val- 
ues of branching ratios obtained in experimental analyzes 
are model dependent. In phcnomcnological studies of B 
meson decays, the isobar model is employed, which has 
a large number of arbitrary phases and relative intensity 
parameters. It is important to use a well-constrained 
model to describe the final state interactions, where the 
constraints can come from theory and analyzes of other 
processes in which the same final states are produced. 
One of our aims is to reduce the number of free param- 
eters by using a model which describes many resonances 
in a unitarized way. 

In Ref. [7], particular attention was paid to the B — > 
/o(980)A" channels which have relatively large branching 
ratios for charmless -B-decays. The amplitudes consid- 



ered in [7] , based on the QCD factorization approach, do 
not include hard-spectator and annihilation terms which 
contain phenomcnological parameters. It was found that 
these amplitudes do not reproduce the B ± — ► / (980)A ± 
and B° — > /o(980)A° branching ratios. Therefore, addi- 
tional terms were included in the decay amplitudes in 
the form of long-distance contributions stemming from 
so-called charming penguins [8]. At the hadronic level, 
these contributions could be associated with intermedi- 
ate DgD^*} states and be understood as being part of 
the intricate final-state interactions in B decays. One 
can expect a sizable contribution from such processes as 
B — > D*D(*) branching fractions are large. The addition 
of charming penguins and of 7T7t and KK coupled-channel 
final state interactions in Ref. [7] allowed to obtain a good 
agreement with the measured B — > /o(980)A branching 
ratios and to reproduce the effective 7T7t mass distribution 
of the BaBar and Belle data. 

Long distance contributions, similar to the charming 
penguin effects, have been considered by Cheng, Chua 
and Soni [9]. Cheng, Chua and Yang have recently 
studied the nature of the scalar mesons Jo (980) and 
a (980) in the decays B -► / (980)A, B -> a (980)7r and 
B — ► ao(980)A using QCD factorization for the weak 
decay amplitudes [10]. Their calculations are for a qq 
state of the /o(980), but implications of a q 2 q 2 picture of 
the /o(980) are also discussed. In their work [10], how- 
ever, they neither consider hadronic long-distance contri- 
butions nor pionic and kaonic final-state interactions. 

In the literature, one finds different results on the rela- 
tive size of the weak annihilation or hard-spectator con- 
tributions to the -B-decay amplitude into two mesons [11- 
13]. A recent study in soft-collinear effective theory [11] 
concludes that the annihilation and chirally enhanced 
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annihilation contributions in charmless B — > M1M2 de- 
cays are real to leading order in Aqcd /wi> (here M\ and 
M2 are non-isosinglet mesons, Aqcd is the QCD scale 
and rrib the 6-quark mass). These contributions consti- 
tute a small fraction of the experimentally determined 
total penguin amplitudes in the case P° — ► K~n + or 
B~ — > K~K° decays and in the present work they are 
not considered. Even upon inclusion of the relatively 
small annihilation and hard-spectator contributions, the 
current QCD factorization results underestimate the av- 
erage branching ratio for B — > pK by a factor of two 
and larger [14, 15] unless certain model parameters are 
strongly modified (see scenarios S2 to S4 in Ref. [15]). To 
be more quantitative, let us cite the branching ratio of 
the B~ — > K~p° decay calculated by Leitner, Guo and 
Thomas [14]. Including fully these two contributions, it 
is equal to 1.54 x 10~ 6 . This should be compared to the 
experimental branching ratio for this decay which lies be- 
tween 3.9 x 10~ 6 and 5.1 x 10~ 6 with a typical error of 
0.5 x 1CT 6 [1, 4]. 

In this paper, we supplement the B — > p(770)°K, 
p(770)° — > 7r + 7r~ decay amplitudes to the B — > 
/o(980)if, /o(980) -» ir+ir- amplitudes derived in 
Rcf. [7]. The experimental branching ratio for the de- 
cay B — ► p(770)°K is of the same order of magnitude as 
for the B — > /o(980)-fT channel. A characteristic fea- 
ture of the B ± — > p(77Q)°K ± decay, however, is the 
very large value of the direct CP violating asymmetry 
of about 0.3 measured by Belle [1] and BaBar [4]. The 
addition of the B — > p(770)°AT decay amplitudes enables 
us to study possible interferences between the S- and P- 
waves in B — > mK decays. These effects can be observed 
in particular in the time dependent CP asymmetries of 
neutral B decays into p(770fK% and / (980)ATg chan- 
nels. In b — > s transitions, new physics contributions 
to these asymmetries have been put forward [16, 17]. 
Nonetheless, one should also take into account the in- 
fluence of meson-meson final-state interactions on these 
observables, among others on the deviation of the asym- 
metry parameter S from the Standard Model expectation 
value sin 2/3. 

In Section II, after a concise summary of the 5-wave 
contributions B — > fo(980)K of Rcf. [7], we derive the 
B — > p(770)°K decay amplitudes from the QCD fac- 
torization approach of Ref. [15]. We introduce a Brcit- 
Wigner vertex function r pn -n(m W7V ) in the B — > p(770)° K 
decay amplitude to account for the tttt final-state inter- 
actions in the P-wave. We complement this amplitude 
with the S- wave contribution B — > f n (980)K obtained in 
Ref. [7] . In Section III, we give our model expressions for 
the observables measured by the Belle and BaBar Col- 
laborations. We briefly discuss our fitting method in Sec- 
tion IV. In Section V our results are compared to the ex- 
perimental 7T7T mass and helicity-angle distributions and 
we give our branching ratio values for B — > p(770)°K and 
B — > fo{980)°K decays. We analyze the S- and P-wave 
interference effects as well as the influence of charming 
penguins on the CP violating parameters 5(m W7r ) and 



A(m,T r7r ) of the time-dependent asymmetry in neutral B° 
decays. Furthermore, we calculate the asymmetry pa- 
rameters S and A and compare our results with experi- 
mental data. In Section VI, we summarize and propose 
further improvements of this work. Finally, for complete- 
ness, we give, in the Appendix the expressions of the 
B~ -» / (980)A- and B~ -» (TT+TT-) S - wave K- decay 
amplitudes. This will also allow us to make a comparison 
with other approaches like that of Ref. [10]. 

II. AMPLITUDES FOR THE B -f tttzK DECAYS 



A. A brief recall of the B — > f (980)K amplitudes 

In Ref. [7], the B decays into ttttK and KKK were 
studied for final-state (7T7r),s and (KK)s pairs interact- 
ing in an isospin zero S*-wave from tttt threshold to about 
1.2 GeV. The QCD factorization approach is the frame- 
work used to describe weak decays of B mesons into a 
quasi two-body state fo(980)K. The effective coefficients 
ai{p){i = 1,...,6) at the renormalization scale p, = 2.1 
GeV were taken from Ref. [18]. Corrections arising from 
annihilation topologies and hard gluon scattering with 
the spectator quark were not included. 

The two-pion and two-kaon rescattering effects arc 
described by the tttt and KK unitary coupled-channel 
model of Ref. [19], which allows to describe the lowest 
scalar- isoscalar resonances /o(600) and / (980) with a 
single matrix of amplitudes. The meson-meson ampli- 
tudes are incorporated into four scalar form factors re- 
sponsible for the production of the 7T7t and KK pairs 
in an isospin zero S- wave from the current uu, dd and 
ss antiquark-quark states. These scalar form factors are 
constrained by the chiral dynamics of low-energy meson- 
meson interactions [20] and provide the link between the 
weak decay amplitudes [18] and the hadronic rescattering 
model [19]. 

As mentioned in the introduction, charming penguin 
amplitudes were included. In a nutshell, these ampli- 
tudes are low-momentum cc loop contributions which are 
formally suppressed by powers of Aqcd /mb though nu- 
merically enhanced by large Wilson coefficients [8]. No 
explicit calculations being available so far, one interprets 
these amplitudes as long distance contributions parame- 
terized by two complex parameters. 

B. The B —>■ pK,p — > (tttt)p decay amplitudes 

The derivation of the B — > /?(770) K weak decay am- 
plitude follows the QCD factorization approach of Bcnckc 
and Neubert [15]. By (7r7r)p we denote isovector tttt 
states in a P-wave and in the following we will write inter- 
changeably p(770)° or (mr)p depending on the context. 
The possible quark line diagrams for negative P-meson 
decays are shown in Fig. 1, where the ir + ir~ final state 
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interaction is graphically schematized by dashed lines ap- 
pended to either a uu or dd state in the diagrams. 

The B~ — > (ir + Tr~)p K~ decay amplitude can be writ- 
ten as 

((7T + 7r~) P K~\H\B~) =Pk-(Ptt- ~ Pw+) T p7rn (m n „), 

where pk is four-momentum of K and p n - and p n + are 
the four-momenta of the negative and positive pions, re- 
spectively. In the 7r + 7r _ rest frame, the decay amplitude 
reduces to 

((■k + tt-) p K~\H\B-) 

= 2A~ r p7r7T (m n „) Ip^UpkI cos6», (2) 



M'x are the mod- 



where \p„\ and |pk| = \J E 2 K {m^) 
uli of the pion and kaon momenta with E^m^^) 
|(M| - - M^/m^n, M B , m„„ and M K being 
the B-meson, effective tttt and kaon masses, respectively. 
Moreover, 8 is the helicity angle between the direction of 
flight of the 7T~ in the tt + tt~ rest frame and the direction 
of the 7T + 7r - system in the B rest frame. 
The function A~ is given by 



A~ 



Gp m p 



1k A 

■ K H 



C p ) 



(3) 



where Gp is the Fermi coupling constant, fx = 
159.8 MeV is the kaon decay constant and the p de- 
cay constant is given by f p = 209 MeV. For the tran- 



sition form factors B 



A^»{M_ 



p and B 



K, we employ 
0.37 and Ff^ n (m 2 p ) = 0.34, respectively. 
The above values for decay constants and transition form 
factors are taken from Table I of Ref. [15]. The functions 
U~ and W~ are defined as 



= A„[ai + a£ - a£ - r (ag -a%+ - a%) 

+ a\ -a c w ] + X t [-al + r{al + a%)-a c w ] (4) 
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The chiral factor is given by r — 2Mg/[(rnt, + m u )(m s + 
m u )], with rrib — 4.46 GeV and where m u = 5 MeV and 
m s = 110 MeV are the u- and s-quark masses, respec- 
tively. The coefficients ai(/i = rrib/2) are taken to be 
the a hI (nK) of Table III in Ref. [21]. In Eq. (4), we 
did not include the annihilation contribution considered 
in Ref. [15]. The charming penguin contribution C p in 
Eq. (3) is parameterized as 



C 1 



(6) 



with P u and Pt being complex parameters [8]. In Eqs. (4), 
(5) and (6), we have used the unitarity condition A c + 



A,, + A* 



0, where At, 



KibKfc: At 



V tb V t * and A c 



FIG. 1: Quark line diagrams for the B~ decay: (a) and (b) 
are tree diagrams, (c) is the penguin diagram and (d) and 
(e) are electroweak penguin diagrams. The spring-like lines 
represent gluon exchange, wavy lines photon or W~ exchange 
and the dashed ones the nn isospin 1 P-wave pair. 



V c bV* s are products of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix elements, to express the amplitudes in 
terms of A u and At . The parameters of the CKM matrix 
are taken from Ref. [22] . 

The corresponding A + function of the B + — > 
(tt~ tt + )p K + decay amplitude is obtained from replac- 
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ing the A u and At values by their complex conjugates A* 
and Aj in Eqs. (4) to (6) and changing the overall sign in 
Eq. (2), 

({ir-7r+)pK+\H\B + ) = 

= 2A+r p7r7r (m 7r7r ) Ip^HpkIcosA, (7) 



where we define 



A+ = -A-(\* u ,\$). 



(8) 



The minus sign in Eq. (8) follows from the CP symmetry 
of the final state. 

Coming to neutral B decays, the B° — > p(770)°-K"°, 
p° — > (-K + n~) p decay amplitude reads 

((n + n-)pK°\H\B } 

= 2A a T p7I ^(m 7T7[ ) IPxIIpatI cos6», (9) 



with 



A = Gf ni„ 



f K A^ p (M 2 K )(U° + C p ) 



+ ,f p F^ K (nip) W u 



(10) 



The different sign in front of the charming penguin term 
C p is due to the dd quark content of the p(770)° in neu- 
tral B decays in contrast to Eq. (3), where in charged B 
decays the uu configuration comes into play. The func- 
tion U° is given by 

U° = X u {-at + al + r[a%-al- « - a$)/2] 

+« - a c 10 )/2} + \ t {at -r(a c 6 - a c s /2) - a\ /2] (11) 

and W° = W~ . Here again we have left out annihilation 
contributions as considered in Ref. [15]. The replacement 
of X u and At by their complex conjugate values in Eqs. (5) 
and (11) leads to the corresponding function A for the 
B° — ► (-k~it + )p K° decays which reads 

((Tr-Tr+) P K a \H\B°) 

= 2A a T p7I7T (m 7T7I ) \p„\\pk\ cos6», (12) 

where one has again to account for an overall sign change 
in the amplitude 



A° = -A°(\* u ,\* t ). 



(13) 



The p — > 7T7T vertex function is chosen to be of Breit- 
Wigncr form, 



C. Complete B — > -kit K decay amplitudes 

We add the S-wave amplitude 05 of Ref. [7] for the 
B — ► (tttt) s K decays to obtain the complete decay am- 
plitude as 



M =a s + a P Ip^ | \p K | cos 9, 



(15) 



where as is given by either Eq. (1) or Eq. (7) of Ref. [7] 
and ap can be read either from Eq. (2) or Eq. (9) for B~ 
decays or for B° decays, respectively. The expression of 
the B~ — ► (Tr + Tr~)sK~ amplitude can also be found in 
the Appendix A (Eq. (A12)). It follows from Ref. [7] and 
from Eqs. (8) and (13) that the CP-conjugate S- and P- 
wave amplitudes as and a p are related to the amplitudes 
as and ap by 

d s = a s {X* u ,\* t ), dp = -a P (X* u , X* t ). (16) 

Thus, the conjugate of M is written as 



M =a s + aplp^WpKl cos6. 



(17) 



In the amplitude as of the present work, the charming 
penguin contribution C (m) of Ref. [7] is replaced by 

C s (m) = -(M 2 B -m 2 )f K F B ^ )s (M 2 K )(\ u S u + \ t S t ). 

(18) 

In Eq. (18), S u and St are two complex parameters and 



J (M|.) = 0.46 [7] is the B-transition form fac- 
tor to a pair of pions in an S-state. Furthermore, the 
constant x, characteristic of the /o(980) — > tttt decay, is 
fixed to 29.8 GeV -1 in agreement with the estimation 
made in Ref. [7]. Here, for the S'-wave, we also use the 
a.ij{p = TOb/2) from Table III of Ref. [21] as effective 
coefficients dj. 

We can substitute the coefficients ai(PP) of the 
.B-decay into two pseudoscalar mesons PP in place 
of the coefficients a,ij(PS) corresponding to the final 
pscudoscalar-scalar PS state for the following reasons. 
As can be seen from Eq. (35) of Ref. [15], the vertex 
corrections, which depend on the scalar light-cone distri- 
bution amplitudes, affect only the aej(PS) coefficients 
of the amplitude for which the emitted meson without 
the quark spectator is a scalar. In our S'-wave ampli- 
tude it is only found in the a^ of Eq. (A5). However, 
aaj(PS) is practically equal to a§j{PP) (see Ref. [10] 
and Ref. [11] cited in Ref. [7]). Any further corrections 
to the a.i are due to hard-gluon scattering contributions 
(a-ijiin) of Ref. [21]) as well as annihilation terms; both 
processes are, however, not included here. 



r p7TTT (^"7T7r) 



9p7T7T 



mt 



- m 2 p + iTpnip 1 



(14) 



where m p — 775.8 MeV, the decay width is T p = 
150.3 MeV and the coupling constant is taken to be 
9 2 p7T j47T = 3m 2 p T p /2pl 



III. DEFINITIONS OF OBSERVABLES 

The B — > nirK amplitude M in Eq. (15) depends 
on the effective mass m W7r and cos 9 which is equivalent 
to the other effective mass variable m V K on the Dalitz 
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plot. The double differential B — ► irirK decay distribu- 
tion then reads 



d 2 T 



d cos 6 dm nir 8 (2tt) 3 Mf 



l-MI 1 



(19) 



Integrating over COS6 1 one obtains the differential B 
ttttK decay distribution 



dT m vw \p w \\p K \ , 

\as\ 



Vir\-\VK\-\a P \ 2 ) (20) 



dm^ 4 (2tt) 3 M| 
and the differential branching ratio is given by 
dB 1 dT 



dm,, 



Tb dm-n 



(21) 



where Ts is the appropriate total width for B + or B°. 

The integration on the Dalitz plot over the kinemat- 
ically allowed m V7r range yields the hclicity angle B — > 
ttttK distribution 



dT 



= A + B cos9 + C cos 2 9. 



dcosd 

The functions A, B and C are defined as 
to^IpttIIpkI 



A 
B 
C 



J m 



\as\ dm-xir. 



(22) 



(23) 



8 (2tt) 3 M| 

8(2tt) 3 M 3 Rc ( a s a p) rfm ^. ( 24 ) 



m^lp^PlpKl 3 

8 (2tt) 3 Ml 



ap\ 2 drrinn- 



(25) 



The CP violating asymmetry for charged B decays is 
defined as 



dT{B- -> Tr+TT-is:-) dT(B+ -> tt" 



if 



dm. 



dT(B~ -» 7r+7r -jC-) aT(£+ 



-7T-if+) 



dm, 



din -j. 



(26) 

Furthermore, in neutral -B meson decays into final CP 
eigenstates / the time-dependent asymmetries are given 

by 



r(£°(f)^/)-r(B°(t)^/) 
r(B0(t)^/) + r(B"(t)^/) 

S sin( Ami) + .A cos( Ami) . (27) 



The mass difference between the two neutral -B eigen- 
states is denoted by Am, 5 is a measure for the mixing 
induced CP asymmetry, and A is the direct CP violating 
asymmetry. The two parameters of the time dependent 
asymmetry, S and A, are 



5 = ,:l .1:1mA,. 



A = - 



1 + IA/l 2 
1-IA/l 2 



|A f l 2 ' 



(28a) 
(28b) 



with 



X f = e 



- -np 



M(B°^f) 
M(B° — > /) 



(29) 



and (3 being the CKM matrix angle. Inserting the cor- 
responding amplitudes M(B° — > /) from Eq. (15) and 
M(B° -» /) from Eq. (17) into Eqs. (28a) and (28b), 
one obtains after integration over cos 9 the integrated S 
asymmetry 



S(m™) = (1-4) x 



x 2Im^e 



-2i/3 



a s«5 + glPTr^lPxl 2 apap 



(30) 



and the A asymmetry integrated over cos 9 is 
A(m nv ) 



as| 2 -|as| 2 + i|p.| 2 | P ^| 2 (|ap| 2 



a P \ 



OS 



±| P .| 2 | P k| 2 (M 2 



IopI 2 ) 
(31) 



IV. FITTING PROCEDURE 



In our fitting procedure we include the experimental 
data of Belle [1, 3, 23, 24] and BaBar [4, 5, 25, 26] Col- 
laborations. They include altogether 18 values of branch- 
ing fractions, direct CP violating asymmetries for the 
charged B decays and the time dependent CP asym- 
metry parameters S and A determined for the B° — ► 
p{77$)K° s and B° -» / (980)ifg decays. These data con- 
stitute the first part of the observable ensemble we fit. 
If one denotes by Z° Kp and Zf 1 the experimental and 
theoretical value of the above mentioned observables and 
by AZ° Kp their empirical errors, then the corresponding 
first part of the x 2 function reads: 



Xi 2 



Wi 



18 

E 

»=i 



zT v - z 



ih 



AZ oxp 



(32) 



Here, W\ = 10 is a weight factor chosen so as to make a 
good balance between this and the second component Xu 
of the total x 2 function. The \\\ function is defined sim- 
ilarly to Eq. (32) but for the observables Z^ being the rnr 
effective-mass and the helicity-angle distributions. We 
set Wn = 1. The experimental distributions have been 
background subtracted according to the figures of the 
Belle and BaBar papers cited above. For completeness, 
we give a list of figures showing the data we used in the 
evaluation of figures 5b, 6c, 6d and 7a till 7f from 
Ref. [23], figures 5c, 6c, 6d from Rcf. [3], figure 3b from 
Ref. [5], and figure 3 from Ref. [4]. The experimental 
data are not absolutely normalized, therefore we have 
adopted an adequate method, described below, to com- 
pare the theoretical differential distributions of Eqs. (20) 
and (22) with the experimental m vv and cos# distribu- 
tions. 
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TABLE I: Charming penguin parameters in the S-wave (S u 
and St) and in the P-wave (P u and Pt). 



Parameter 


Value 


\S U \ 


0.15 ± 0.10 


Arg b u 


1.90 ± 0.71 


\St\ 


0.020 ± 0.002 


Arg St 


-0.26 ±0.21 


1^1 


1.09 ±0.21 


Arg P n 


-0.98 ±0.12 


\Pt\ 


0.065 ± 0.002 


ArgPt 


-1.56 ±0.08 



For both differential distributions in x — m V7I and x = 
cos 9, we define the theoretical number of events in a bin 
Xi by the integral 



Y th ( Xi ) =Afr 



dTjx) 
dx 



dx. 



(33) 



The differential model distribution dT(x)/dx is just that 
given by either Eq. (20) or Eq. (22) and the integration 
range of the kinematic variables is the bin width [xn , 2^2] • 
Then the normalization coefficient TV of a given experi- 
mental distribution y exp is defined as a sum over all cho- 
sen bins: 



TV 



12i=l Y C xp(Xj 

~x 2 dT{x) 
dx 



(34) 



dx 



The sum of experimental events Y exp (xi) over Xi and 
the integration of dT(x)/dx over x is done in the an- 
alyzed range [Xi, X2] of interest, which must be within 
the kinematically allowed range for the variables and 
cos 9. Since we concentrate ourselves on the 7T7t effective 
mass range where the p(770)° and /o(980) resonances 
contribute the most, we take in our fits X\ = 0.60 GeV 
and X 2 = 1.06 GeV. 

The experimental branching fractions for charged and 
neutral B -> p(770)°K and B -> ,f (980)K decays arc 
obtained with the isobar model, in which the -B-decay 
amplitude is a sum over the contributions of different 
two-body resonances. These are usually parameterized 
in terms of Breit-Wigner functions with the exception of 
the /o(980), for which the Flatte formula is used. Af- 
ter fitting the intensities of the resonances, a particular 
branching fraction for a given resonance is calculated as 
an integral over the fully available phase space on the 
Dalitz plot. In order to compare these experimental 
branching fractions with those for limited effective mass 
ranges, we have calculated appropriate reduction coeffi- 
cients as the ratio of the integral over the limited phase 
space to that over the fully allowed one using the phe- 
nomenological amplitudes of the experimental analyzes. 
They are equal to 0.71 in the case of the p(770)° and a 
m W7r range 0.66 - 0.90 GeV and to 0.69 for the / (980) 



which dominates in the 0.90 — 1.06 GeV range. In fit- 
ting our model to the data, these coefficients reduce the 
branching fractions and their errors by the same factor. 

Experimental cuts on the effective pion-kaon mass in- 
fluence the 7T7T effective-mass and the helicity angle distri- 
butions. A typical veto cut is experimentally used around 
the Kit effective mass, close to the D-meson mass. In 
our calculations, we have applied the cuts specified by 
the Belle and BaBar Collaborations in their analyzes of 
the B decays into 7r7r.fr. 

Our model has four complex parameters: P u and P t 
in the P-wave amplitude and S u and St in the S'-wave 
amplitude. With these parameters we shall describe more 
than two hundred data values. 



V. RESULTS 

We have performed a global fit to the available data 
measured by BaBar and Belle. These consist of 204 data 
points describing the 7T7t mass and angular distributions 
and of 18 observables enumerated in Section IV. The \ 2 
values are the following: Xj/Wi = 9.8 and \ii = 336.3. 
The values of the eight real parameters are listed in Ta- 
ble I along with the parabolic errors from the MINUIT 
minimization procedure [27]. 

In Table II, the values of branching fractions, direct 
and time-dependent CP asymmetries are presented along 
with the corresponding experimental ones. The calcu- 
lated observables for a given channel are integrated over 
the range indicated in this table. Both the BaBar 
and Belle values are in agreement within their error bars 
and are well reproduced in our model. The theoretical 
errors stem from the parameter errors given in Table I. 

In Subsections VA and VB, we present the 7T7t mass 
and helicity-angle distributions for charged and neutral 
B — > ir + ir~K decays, respectively. Interference effects 
are studied in Subsection VC and some discussion on 
the results is given in Subsection VD. 



A. The 7T7T mass and helicity-angle distributions in 
— > tt + tt~K decays 

In Fig. 2 the 7T7t effective-mass distributions of our 
model are compared to the Belle data [23]. In the fit 
to these distributions, the background corrected data 
have been used. One sees from Figs. 2a and 2b that our 
model describes rather well the 7r + 7r~ spectra measured 
by the Belle Collaboration separately for the B + and 
B~ decays. Our theoretical curve depicts, as in previous 
work [7], a prominent maximum near 1 GeV, to which 
now adds a less pronounced peak at about 770—780 MeV. 
Remarkably, one observes a clearly visible asymmetry 
in number of events between the B~ — ► K~-k + it~ and 
B+ -> K+ir+ir- decays for the ^(770)° and / (980) re- 
gions. At lower masses about 500 MeV, our model 
also produces a broad maximum which we attribute to 
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TABLE II: Average branching fractions B in units of 10 , asymmetries Acp, and the asymmetry parameters S and A of our 
model compared to the values of Belle and BaBar Collaborations. The experimental branching ratios and their errors have 
been multiplied by 0.71 for B — + pK channels and by 0.69 for B — > foK channels (for explanation see Section IV). 

Observable Channel ir + ir~ mass Our model Belle BaBar 
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0.02 ± 0.02 


-0.08 ± 0.08 
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B 


p°K° 


(0.66,0.90) 


3.49 ± 0.20 


4.35 ± 1.05 


[3] 


3.62 ± 1.14 


[5] 


B 


foK 


(0.90, 1.06) 


3.86 ±0.14 


5.24 ± 1.28 


[3] 


3.79 ±0.62 


[5] 


S 


TT + %-K^ 


(0.66, 0.90) 


0.11 ± 0.11 






0.17 ±0.58 


[26a] 


A 


TV + TT-K^ 


(0.66, 0.90) 


0.01 ±0.06 






-0.64 ±0.48 


[26a] 


S 


■K + -R-K% 


(0.89, 1.088) 


-0.67 ±0.05 


-0.47 ±0.37 


[24a] 






S 


■K + TT-K% 


(0.86, 1.10) 


-0.63 ±0.05 






-0.95 ±0.29 


[25] 


A 


■K + TT-K% 


(0.89, 1.088) 


-0.11 ±0.07 


-0.23 ± 0.26 


[24a] 






A 


ir+n-K° s 


(0.86, 1.10) 


-0.11 ±0.07 






0.24 ± 0.34 


[25] 



the a or /o(600). By the same token, we note that in fits 
to their data, the Belle Collaboration has not included 
the a resonance which is thus buried in the non-resonant 
background. Let us remark that this part of the m„ 
spectrum has not been fitted in our calculations. As ex- 
plained in the previous section IV, we have intentionally 
chosen the lower mass limit to be equal to 0.6 GeV. 

In Fig. 3, we compare the results of our model for the 
cos 9 distribution in the vicinity of the p(770)° (Fig. 3a) 
and of the / (980) (Fig. 3b) with the Belle data [23]. The 
experimental cos 9 distribution is fairly well reproduced 
given the fluctuation in events of the data points. The 
behavior of the angular distribution such as seen near 
the /o(980) mass is typical for a substantial interference 
pattern between the S- and P-wave. The former gives 
a constant contribution while the latter is a symmetric 
cos 2 9 distribution observed in the p(770)° mass range 
(see Eq. (22)). The S — P interference term, which is 
proportional to cos 9, is clearly seen in Fig. 3b. This 
demonstrates the presence of the P-wave contribution in 
the / (980) region. 

B. The 7T7T mass and helicity-angle distributions in 
B° -> tt+tt-K decays 

The m W7r effective mass distribution for neutral B de- 
cays is plotted in Fig. 4 along with the results obtained by 
BaBar [5] . Within the error bars our model also describes 
this distribution well, which bears the same features 
concerning the p(770)° and /o(980) as the charged B de- 
cay m TTr spectrum. In Fig. 5a and 5b the cos 9 distribu- 
tions in the proximity of the p(770)° and /o(980) are com- 
pared with the experimental data from Belle [3]. While 
the p(770)° range is dominated by the P-wave contribu- 
tion, in the /o(980) range one observes once again the 
interference pattern. The sign of the interference term 



is, however, opposite to that apparent in Fig. 2b. It is 
well explained by our model in which the P-wave and S- 
wave amplitudes for both charged and neutral B decays 
are closely related (see Eqs. (3) and (10) where the signs 
of the charming penguin terms are reversed). 



C. Interference effects 

In this subsection we study the interplay between 
the S- and P-waves of the decay amplitudes. Fig. 6 
shows the angular dependence of the two CP violating 
asymmetries S and A calculated at the effective masses 
m^t = 775.8 MeV and m„ = 980 MeV using Eqs. (28a) 
and (28b) and (29). These numbers are chosen to demon- 
strate the behavior of the asymmetries about the p(770)° 
and /o(980) resonances. As seen in Fig. 6a and 6b, the 
functions S and A are not constant. This is a clear 
demonstration of the interference between S- and P- 
waves. The sharp changes in behavior of the asymme- 
try S at the p(770)° mass are shown in Fig. 6c, where 
the straight dashed and dotted lines are the S values for 
separate S- and P-wave contributions, respectively. At 
cos 9 = 0, the P-wave amplitude vanishes, therefore the 
strongest variations of the asymmetry are observed near 
9 = it/2. A similar picture for the asymmetry A is shown 
in Fig. 6d. 

The asymmetries S and A at = 980 MeV behave 
smoothly near cos 9 = since there the S- wave domi- 
nates. A decrease of A as well as an increase of S as 
functions of cos 9 are due to the P-wave amplitude com- 
ponent, which is still non-negligible in the /o(980) mass 
range. 

Recall that the helicity angle dependence of the asym- 
metries can be transformed into the functional depen- 
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m^CGeV/c 2 ) (GeV/c 2 ) 

FIG. 2: The 7r + 7r~ effective-mass distributions (a) in B~ — > 7r + 7r~if ~ and (b) in B + — > 7r + 7r _ .K' + decays. The data are taken 
from the Belle Collaboration [23]. The solid lines represent the results of our model. 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 

cos 9 cos 9 

FIG. 3: Helicity-angle distributions in B — ► -k + tx' -K* (a) in the p(770) mass region (0.6 GeV < m nn < 0.9 GeV) and (b) 
in the /o(980) region (0.9 GeV < m^^ < 1.06 GeV). The data points are from Ref. [23] and the solid lines denote our model. 
The S-wave contribution is plotted with dashed lines, the P-wave with dotted lines and the interference term with dot-dashed 
lines. 



dcncc on the tt K effective mass 



m v - K = yml + M\ + m ww E K (m nw ) + 2 |p w ||pjf| cos0, 

(35) 

where is the pion mass. Thus, the above interference 
effects can also be experimentally studied by examination 
of particular regions of the Dalitz plot. 

In Fig. 7, the effective irir mass dependence of 
the asymmetries integrated over cos 8 are plotted (sec 
Eqs. (30) and (31)). The most pronounced effect is an 
intermediate change of sign of S near the p(770)° res- 



onance. We stress that if the S- or P-amplitudes were 
dominated by only one weak amplitude proportional to 
At then the asymmetry S shown in Fig. 7a would sud- 
denly jump from — sin 2/3 w —0.7 at the 7T7t threshold 
to + sin 2/3 w +0.7 near the p(770)° mass and abruptly 
drop down to — sin2/3 w —0.7 at the /o(980) mass, as 
can be inferred from Eq. (30). A smooth behavior of S is 
explained by the mass dependence of the S'-wave scalar 
form factors and by the finite width of the p(770)°. As 
seen in Fig. 7a, our model is in agreement with the data. 
Remarkably, there is a particularly large departure from 



9 




m (GeV/c ) 

7131 v ' 



FIG. 4: The tt + tv~ effective-mass distributions in B { — > 
■k + it~K% decays. Data points are taken from the BaBar 
Collaboration [5] while the continuous line results from our 
model. 



the + sin 2(3 value near the p(770)° mass. This can be 
explained by a substantial value of the weak part of the 
amplitude in Eqs. (6) and (10) proportional to X u . In 
this context, bear in mind that the P-wave penguin pa- 
rameter \P U \ is much larger than the S*-wave parameter 
\S U \ (see Table I). 

The functional dependence of the asymmetry A on 
m vv is depicted in Fig. 7b. It is not as strong as for 
the S asymmetry, but certainly A is not equal to zero. 
In the /o(980) range, A is close to —0.1. Hence, it is 
comparable to the values of A found by Belle and Babar 
in B° — > K + ir~ decays [28, 29]. In the former case, how- 
ever, the experimental errors are yet too large to claim 
the non-zero value of A (see Table II). 

We conclude this section by remarking that in pen- 
guin dominated decays, such as B — > pK,u)K..., tree di- 
agrams are either CKM suppressed or CKM and color 
suppressed. Still, we observe that for p(770)°K final 
states an accidental cancellation occurs in the A„ terms 
of Eqs. (4) and (11) due to the combinations a\ — ra§ and 



— rag , where r 



1. 



« a 4 ' . A similar cancellation 
takes equally place in the B — > /o(980)-ftT channel [7]. 
Therefore, the 'tree pollution' is considerable and one 
may, at this point already, expect deviations from the 
value S s=s sin 2(3 associated with the pure tree diagram 
for the B° — > J/ipKs or with the penguin diagram for the 
B° — ► c/)K decays. The tree diagram contribution to the 
p(770)°K state, however, is not sufficient to explain the 
important departure of S from the above value. It is the 
aforementioned long-distance contribution \P U \ to the A„ 
term combined with interference effects that determine 
the functional behavior of the S asymmetry parameter. 



D. Discussion of the results 



We stress the importance of charming penguin terms. 
Without them it is not possible to obtain a good agree- 
ment of the theoretical branching ratios with experimen- 
tal data. If all the penguin parameters are set equal to 
zero then the model branching ratio is underestimated by 
a factor of about 5 for the B° — > fo(9SO)K° decay and 
by 3.5 for the B° -> p(770)°K° decay. For the charged B 
decays the theoretical branching ratios are too small by 
factors of about 2.4 and 7 for the f (980)K and p{770fK 
channels, respectively. 

It has, however, been found in Rcf. [10] that the 
branching ratios for B — > fo(980)K decays can be repro- 
duced in the QCD factorization approach provided that 
one uses a large scalar decay constant fj o of 370 MeV. 
The work of Ref. [10] is based on a two-body B-decay am- 
plitude. Our approach considers a three-body B-decay 
amplitude by taking into account the /o(980) disintegra- 
tion into two pions or two kaons and also by introducing 
final-state interactions between them. A relation between 
the two- and three-body £>-decay amplitudes, elucidated 
in the Appendix A, leads in our case to a value of 94 MeV 
for fj o as seen in Eq. (A14). The difference with the 
value of Ref. [10] by a factor of about 4 explains partly 
why we underestimate the B — > fo(980)K branching ra- 
tio if long-distance charming penguin contributions are 
absent. 

We have studied the sensitivity of our results on the 
scalar form factors r^l* 2 { m -n-K) which depend on the not 
very well known low-energy constants of the chiral per- 
turbation theory. The possible variation of their values 
(see Refs. [20] and [30]) can lead up to a multiplicative 
factor of about 1.25 for r™!* 2 (m T7r ) in the /o(980) range. 
Note however, that the sensitivity to such changes is lim- 
ited due to the constant value of the product x\^i ( m / )l 
given by Eq. (A15). If, for example, r™(m/ ) is multi- 
plied by a factor of 1.25, the \ value is to be divided by 
1.25 and consequently the value of fj o is also multiplied 
by 1.25 (see Eq. (A14)). Thus, an increase of the scalar 
form factor rj(m/ ) will only enhance the contribution 
of the Q(m nn )V term in Eq. (A12). 



Our B — ► (tttt)sK amplitude depends also on the 
Fq^^ s (Mj^) transition form factor which is a model 
dependent quantity. Decreasing its value from 0.46 (see 
Subsection II C) to 0.25 (value quoted in Ref. [10]) will 
lead to a readjustment of the charming parameters S u 
and St while fitting the experimental data. The qual- 
ity of the fit is good and similar to that presented here. 
We have also obtained a good fit using the values of the 
low energy constants determined in Ref. [30] (the values 
of the first line of Table I therein) for the scalar form 
factors together with F^ ( ™ )s {M 2 K ) = 0.25. Our con- 
clusions are not changed by these possible modifications 
of our input. 
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FIG. 5: As in Fig. 2 but for B° ir + ir~K% decays and for data from Belle [3]. 



VI. CONCLUSIONS AND OUTLOOK 



In our studies of B — > tt + tt~K decays, we have ex- 
tended the work of Ref. [7] by adding the P-wave con- 
tribution to that of the S-wave in the two-pion final- 
state interactions. The total amplitude contains the res- 
onant channels B -> p(770)°K and B -> f (980)K. This 
work goes beyond the usual approach of calculating two- 
body branching ratios of the B decays, as we have taken 
into account parts of the final state interactions between 
particles in three-body decay channels. In this paper, 
we have analyzed the interactions between pairs of pi- 
ons in the Dalitz plot section with the 7T7t effective mass 
range from threshold to about 1.2 GeV. In the 5- wave 
the rescattering or the transition amplitudes between two 
pions or two kaons are described by the unitary coupled- 
channel model of Ref. [19]. 

There are two components in the weak transition am- 
plitudes. The first term is derived within the factoriza- 
tion approximation with some QCD corrections without 
hard-scattering and annihilation terms. The second con- 
tribution, called charming penguins, is a long-distance 
amplitude originating from penguin-type diagrams with 
c-quark or it-quark loops. Four complex charming pen- 
guin parameters, common for B + , B~ , B° and B° decays, 
have been introduced and fitted to numerous experimen- 
tal data, including the 7T7T effective mass and helicity an- 
gle distributions, branching fractions, direct asymmetry 
Acp and time dependent CP violating asymmetry pa- 
rameters S and A. Our theoretical model reproduces 
well the experimental results. Without the charming 
penguin amplitudes, the model branching fractions of 
B -> f (980)K and B -> ^(770)°^ decays are several 
times smaller than the experimental values. In the m nn 
range below 1.2 GeV, the two main resonances p(770)° 



and /o(980) give rise to important interference effects, 
best visible in the helicity-angle distributions. Thus, in 
Figs. 3 and 5, we notice sizable interference terms of op- 
posite signs comparing the B ± — ► tt + tt^K ± decays with 
the B° -> tt+tt-K ones. 

In our model, the direct CP asymmetry in B ± — > 
p^K^ decays is in agreement with the large experimen- 
tal value measured by the Belle and BaBar groups. The 
parameter S of the time-dependent CP asymmetry for 
B° — ► p°Kg decays is significantly smaller than the value 
sin 2/3 expected in case of a full dominance of the weak 
decay amplitude proportional to At. This is related to a 
large value of the w-penguin parameter P u in the P-wave 
(see Table I). The numbers S u , S t , P u and P t should 
be treated as phcnomcnological parameters which we at- 
tribute to the long-distance penguin contributions. In 
our opinion, however, they can also contain contributions 
from other processes omitted by us like annihilation or 
hard-spectator interaction terms. 

In this approach, all the resonances appear in a natural 
way as poles of the meson-meson amplitudes. No arbi- 
trary phases nor any relative intensity parameters for the 
two discussed resonances p(770)° and /o(980) are needed. 
This is in contrast to the usual phenomcnological ana- 
lyzes in terms of the isobar model applied to fit the Dalitz 
plot density distributions. Let us remark that the scalar 
resonance / (600), as a pole of the 5- wave amplitude, is 
also present in the 7T7T effective mass spectrum. However, 
it has not been included in the phcnomcnological models 
of the Belle and BaBar Collaborations [2, 4, 23]. 

So far, we have restricted our analysis to an effective 
7T7r mass of about 1.2 GeV. In further studies, one can 
extend not only this mass range to include further scalar 
and vector resonances but also treat the other variable 
m^K on the Dalitz plot and the accompanying vector 
resonances K* . Given our findings for the /o(980) and 
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FIG. 6: Helicity-angle dependence of the CP violating asymmetries <S and A calculated for B° — > tt + tv~ K s decays at m„-„- 
equal to m p = 775.8 MeV (solid lines in (a)-(d)) and to m^ ( 980 ) = 980 MeV (dashed lines in (a) and (b)). Separate S-wave 
(dashed lines) and P-wave results (dotted lines) are shown in (c) and (d) for m nn = m p . 



p(770)° resonances, we expect many more interference 
effects at larger effective 7T7t masses and also for the nK 
effective mass range. 
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APPENDIX A: AMPLITUDES FOR THE 

B~ -> / (980)AT- AND B~ -> (ir + n-) s K- 
DECAYS 

We give here the expressions for the two-body 
{f K-\H\B~) and the three-body ((tt+tt-) s K-\H\B J ) 
decay amplitudes and show how they are related to each 
other. In the QCD factorization approach, adding charm- 
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FIG. 7: Effective-mass dependence of the CP violating asymmetries, integrated over cos 9, S in (a) and A in (b) obtained with 
our model (solid lines). The Belle Collaboration results [24b] for the /o(980) resonance are plotted as circles and the BaBar 
data for the / (980) [25] and ,0(770)° [26b] as squares. 



ing penguin terms, the former amplitude can be written 
as 

(f K-\H\B-) = ^ [P(m fo )U + C s (m f0 ) + 
+ Q(m fo )V + C S (M K )] . 

If electrowcak penguin and annihilation terms are not 
included, 

P(m fo ) = f K {M% - m 2 fo )F B ^° (M 2 K ) , (A2) 



Q(m fo )-- 



2{M ~ Ssl %ll-M 2 K )Fr K (rnl) , (A3) 



m b -m s 



U=\ u {a 1 +a 4 i -a c 4 + (a c e -a%)r}+ \ t (a c 6 r-a c 4 ) , (AA) 



V = X U K - <) + *ta c 6 



(A5) 



where m/ is the / (980) mass and F Q B (M|-) and 
F Q B ^ K (m 2 fo ) are the B -> / (980) and B —> K transi- 
tion form factors. The scalar decay constant fj o can be 
defined by the matrix element 

(f \ss\0) = m f J% . (A6) 

The three-body decay amplitude reads 

{(n+n-) s K-\H\B-) 

= ^\fl{ [P(m^)U + C s (m^)] r»_(m™)- 



[Q(m^)V + C S (M K )] ^(m^)}, 



(A7) 



where TJ o7r Jm^ 



and T^ oVV (m W7V ) are the non-strange 



and strange vertex functions, respectively. They describe 
the decay of the /o(980) into two pions. The factor 
\/2/3 is the Clebsch-Gordan coefficient which relates the 
isospin-0 S'-wave |(7T7r)s) to the |(7r + 7r _ )s) one. The ver- 
tex function ^f onn {m nn ) is defined through the relation 
between the matrix elements 



((7T7r)s|s S |0) - T s fown (m wn )(f \ss\0) . (A8) 

Replacement of ss in Eq. (A8) by hn — -^(uu+dd) yields 
the corresponding Vj o7r7r (m n7r ) vertex function. The ma- 
trix element ((7T7t)s|ss|0) is related to the strange scalar 
form factor Tf (m TW ) by [20] 



((7T7r) s \ss\0) = V2B Tl*(m nv ) 



(A9) 



where Bo = — (0\qq\0)/f 2 is proportional to the quark 
condensate, f n being the pion decay constant. Let us 
note, by this occasion, that in Eq. (A9) the charge con- 
jugation appears since the form factor Tf is defined as 
being proportional to the matrix element (0|ss|(7T7t)s). 
One must also correct Eq. (11) of Ref. [7] in which T™' 5 
should be replaced by T"' s *. Using Eqs. (A6, A8, A9) 
one finds 



T %^i m ^) = X T{*(m nv ) , 
where the constant x 1S defined by 

V2B 



X = 



™fof 



(AW) 



(All) 



fa 



We assume rj oJrjr (m 7r7r ) to be related to the non-strange 
scalar for factor r™*(™ 7r7r ) by the same constant x as in 
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Eq. (A10). Here, r™*(m 7rT ) = ((wn) s \nn\0) / V2B simi- 
larly to Eq. (A9). Making the appropriate replacements 
for T s foW7r (m nw ) and T n fo ^{m^) in Eq. (A7), we finally 
arrive at 

{(n+n-) s K-\H\B-) 

= ^\fl{x [P(m^)U + C s {m^)} ry*(m mr )+ 

+ [Q(m nn )V + X C S (M K )] r;*(m™)} , 

(A12) 

where the relation 



Q(m V7r ) = <3(m 7r7r )x 



0413) 



holds. Thus, Eq. (A12) is identical to Eq. (1) of Ref. [7]. 
From Eq. (All) it follows 



V2B 
m foX ' 



(AU) 



With x = 29.8 GeV- 1 , B = ml /2m and m = 5 MeV, 
m/ = 980 MeV one obtains jj a = 94 MeV. In our model 



fj o is determined by the /o(980) properties since [7], 



m /o r tot (/ ) |r«(m /o )| 



(A15) 



here r to t(/o) is the total /o(980) width. The constant \ 
enters in the description of the decay of the scalar reso- 
nances into two pions as seen from Eq. (A10) and com- 
ments below Eq. (All). Furthermore Eq. (A15) shows 
that x has the same role as the coupling constant g p7T7r 
of Eq. (A14). 

As can be seen from Eqs. (Al) and (A6) the two-body 
part of QCD factorization amplitude without charming 
penguins is similar to the corresponding one of Ref. [10] 
when excluding electroweak penguin and annihilation 
contributions. Note, however, that we disagree with the 
relative sign between the P(m T7r )C/ (Eqs. (A2) and (A4)) 
and Q(m T7r )V (Eqs. (A3) and (A5)) terms. Our relative 
sign agrees with that of Ref. [31]. 
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